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ABSTRACT 



The performance characteristics of a residential solar assisted water- to- 
air heat pump heating system are analyzed for State College, Pennsylvania. A 
realistic residence and solar assisted water-to-air heat pump system are modeled 
for this northern climate using the transient simulation computer code TRNSYS 
developed by the University of Wisconsin. The system is studied over a one 
month winter period, December, using actual hourly weather data for State 
College. The system is analyzed for both the cloudiest and clearest December 
weather recorded in the last 30 years. The collector area and storage tank 
capacity are varied and the effects on system performance are studied. Other 
parameters varied and studied are the collector/storage tank heat exchanger 
effectiveness and the room temperature. 
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ABSTRACT 



The performance characteristics of a residential solar 
assisted water-to-air heat pump heating system are analyzed 
for State College, Pennsylvania. A realistic residence and 
solar assisted water-to-air heat pump system are modeled for 
this northern climate using the transient simulation 
computer code TRNSYS developed by the University of 
Wisconsin. The system is studied over a one month winter 
period, December, using actual hourly weather data for State 
College. The system is analyzed for both the cloudiest and 
clearest December weather recorded in the last 30' years. 
The collector area and storage tank capacity are varied and 
the effects on system performance are studied. Other 
parameters varied and studied are the collector/ storage 
tank heat exchanger effectiveness and the room temperature. 
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CHAPTER 1 



1 • 1 In troduc ti on 

With the energy crisis worsening and conventional 
heating fuels becoming scarce and expensive a need exists 
for the development of alternative heating fuela and 
systems* Solar energy which is **free** and readily available 
is a promising source of building heat. Unfortunately for 
northern climates direct solar heating of buildings can 
provide only a small portion of the total heating 

requirements without extremely large collector areas and 
heat storage capacl ty ( 1 ) .* The solar energy can however be 
used to augment conventional heating systems or used to 
improve the performance of the system^ as with a heat pump. 
The heating capacity of the heat pump decreases as the 
evaporator temperature (i.e. the outdoor ambient air 

temperature) decreases. At the same time the heating load 
of the residence increases. The point at which the heat 
pump capacity equals the heating load is called the balance 
point. The heat pump is an economical heat source as long 
as it can operate at a temperature above its balance point. 
In cold northern climates the evaporator inlet temperature 
of the standard air-to-air heat pump installation is often 
below the balance point. For water-to-air heat pumps an 



* Numbers in brackets refer to references listed at the end 

of this paper. 
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existing warm water source (temperature above the balance 
point) for the heat pump evaporator is difficult to find in 
northern climates particularly for residential 
installations. Solar energy can be used easily to provide 
this warm temperature source for the heat pump evaporator 
since only a low grade temperature (50-100 degrees F) is 
r equi red • 

Over the past several years many papers and studies 
have been done on the evaluation of solar energy heat pump 
systems. The earliest research on solar assisted heat pumps 
was begun by the American Gas and Electric Service 
Corporation in 1949 in a 5 year study which used a high 
temperature storage system with the heat pump coupled 
directly to the col le ct or ( 2 ] . In 1952 R. C. Jordan and S. 
L. Threlkeld published a paper in which the potential of 
this system on the basis of some experiments was 
desc r lb ed [ 3 ] . Since that time the solar assisted heat pump 
system has been improved upon and extensively studied. The 
effects of geographical location and variations in the 
system design or component design have been addressed in the 
p re sent literature (4-61. 

The study described in this paper will analyze the 
performance of a solar assisted water-to-air heat pump 
residential heating system in State College, Pennsylvania 
(latitude 40.5 degrees N). A realistic residence and solar 
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assisted water^to^alr heat pump system are modeled to study 
the performance characteristics of such a system In a 
northern climate. The transient simulation computer code 
TRNSYS developed by the University of Wisconsin and as 
modified at the Pennsylvania State University is used. The 
system is studied over a one month winter period, December, 
using actual hourly weather data for State College. The 
system is analyzed for both the cloudiest and clearest 
December weather recorded in the last 30 years. The 
collector area and storage tank capacity are varied and the 
effects on system performance are studied. Other parameters 
varied and studied are the col le c tor /s tor age tank heat 
exchanger effectiveness and the room temperature. 
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CHAPTER 2 



2.1 Description of Residential Structure 

For Investigation and modeling purposes the solar 
assisted heat pump system was assumed to be Installed In a 
single family residence. It Is a two story structure with a 
full basement and attached two car garage. The living area 
Is approximately 2200 square feet with 8 foot ceilings 
throughout. The solar collectors are located on the south 
slope of the roof and are tilted at 50.5 degrees (latitude + 
10 degrees), the recomended [2] optimum angle for heating 
applications. A floor plan of the residence Is shown In 
Figure 1. The residence Is located In State College, 
Pennsylvania, a small university oriented city situated In 
the center of the state. 

An energy conserving design has been utilized to 
minimize the required heating load. The house walls are 
built of wood frame construction with four Inch face brick. 
All windows are double pane glass and a storm door Is 
Installed on the front door. Three and one half Inches of 
R-1 1 Insulation Is Installed In the walls and six Inches of 

r 

R-22 Insulation Is placed In the second floor celling. A 
more complete list of the construction details can be found 
In Appendix A. 
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plan of residence utilizing a 
assisted heat pump system In 
State College, Pa. 



FIGURE 1 



F 1 oo r 
solar 



In order Co size and select Che appropriate heat pump 
It was necessary to calculate a design summer cooling load 
as well as Che winter heating load* The heating load was 
calculated using the general procedures set forth by 
ASHRAE[3] using ASHRAE tabulated data* The resulting design 
heat loss was approximately 32222 Btu/hr* From this the 
building UA (thermal transmittance X exposed area) was found 
to be 511 Btu/hr-F which Indicates a very well insulated 
house. The cooling load was determined using the ASHRAE 
simplified procedures for residential ca Ic u la t Ions ( 4 ] and 
ASHRAE tabulated dat?. The cooling load requirement was 
found to be approximately 1.4 Cons. Summary calculation 
sheets for both the heating and cooling loads can be found 
In Appendix A. 

2.2 System Description 

A schematic diagram of the residential solar assisted 
heat pump heating system Is shown In Figure 2. The major 
components of the system are a flat place solar collector, 
storage tank and a water-to-air heat pump. Other equipment 
includes pumps, heat exchanger, a relief valve and auxiliary 
electric heating which is an integral part of the heat pump. 
The major components are discussed In more detail in the 
following paragraphs. 
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FIGURE 2 - Schematic of solar assisted 
heat pump heating system 



Flat Plate Collector - The flat plate solar collector 
chosen for the simulation was a Pittsburg Plate Glass(PPG) 
Model CD414. The PPG collector features weather sealed 
construction with two tempered glass cover plates and an all 
copper absorber plate with a non-s e lec t 1 ve coating. The 
collector is especially designed for applications where the 
collector area is limited or where large temperature 
differentials exist between the ambient air and the absorber 
plate. Each collector panel is approximately 20 square 
feet. The collector array uses a 40% ethylene glycol-water 
solution as the circulated heat transfer medlumit This 
provides protection from freezing to -12 degrees F. 
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ADDITIONAL PARAMETERS 



Wind Velocity: 0-5 mph 

Fluid: 40% ethylene glycol solution 

Cp - .792 Btu/lb-F 

Flowrate • .025 gpm/sqft 

Tube Diameter • .5 in 

Tube Spacing • 3 in 




FIGURE 3 - Efficiency of the PPG flat plate 

solar collector 



The performance of this flat plate collector can be 
described by a number of parameters. Perhaps the roost basic 
performance parameter is the collector efflciency( ) which 
for the PPG collector is 



-.69-. 7 ( (Tf i-Ta) /HT] ( 1 ) 

where Tfi-inlet fluid t erop e r a t u r e ( d e g r e e s F) 
Ta -ambient air temper ature ( degrees F) 
HT -incident radiation per unit area on 
the collector surface(Rtu/hr-sqft) 
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The efficiency equation is shown graphically in Figure 
3 along with a number of ocher parameters describing the 
construction of the PPG flat plate collector. These 
parameters include the transmittance ('tO of the cover, the 
absorptance (o<) of Che absorber plate and the emmlsalvlty (£) 
of the absorber plate. 

Three other useful parameters describing the PPG 
collector are the Collector Efficiency Factor(FR), the 
Collector Overall Energy Loss Coef f icient (UL) and the 
Collector Geometry Efficiency Factor(F'). These parameters 
are related to the Rate of Useful Energy Collected per Unit 
Area of Collector (qu) according to the following expression 
developed by Hoctel and Woertz[5] and extended by 
Whlllicr{6] and Bliss(7]. 

qu-FRtHT(r^)-UL(Tf 1-Ta)] (2) 

Since «qu/HT equatlon(2) can be rewritten as 

-FR(To< )- [FR*UL(Tf i-Ta)/HT) (3) 

Comparing this to the collector efficiency equation for 
the PPG collector, it follows that 



FR(T«<)-.69 and FR*UL-.7 



(4) 




I 





Substituting known values for the transmittance and 



absor p tance 



FR-, 69/( . 72) (.95)-!. 0 



thus 



UL-.7 Btu/hr-sqf t-F 



(5) 

( 6 ) 



Finally, using UL, the Collector Geometry Factor(F 0 
can be obtained from graphs developed by Duffle and 
Beckman(8]. The value of F' for the PPG collector is .99« 



Storage Tank - An insulated steel underground tank 
containing water is used for thermal storage. The overall 
thermal t ransmit tance (U ) of the tank is .05 Btu/hr-sqf t-F . 
The tank is cylindrical in shape with a constant length to 
diameter ratio of three(as the volume is varied in the 
simulation). A tube-in-shell heat exchanger with an assumed 
effectiveness of 70Z is used to transfer heat from the 
collector loop to the storage tank. 

C. Wat er-t o-Air Heat Pump - An American Air Filter(AAF) 
water-to-air heat pump Model HW19 (1*5 tons) was selected 
for the residence. As is standard practice the heat pump 
was selected to meet the cooling load requirements of the 
building. It is recognized that in harsh weather the heat 
pump will be unable to provide 100% of the heating load 
therefore the AAF heat pump is equipped with a parallel 



auxiliary electric resistance heater to provide the deficit. 
The AA F Model HW19 water-to-air heat pump is rated to 
provide 600 CFM of heated air to the house. The water flow 
on the water side of the heat pump is 5 gpm. Table 1 is the 
heating performance characteristics of the Model HW19 for 
maintaining the room temperature at 70 degrees F. 



TABLE 1 - Heating performance characteristics 





of 


AAF Model HW19 


hea t 


pump 




EVAP 


HEAT 


HEAT 








INLET TEMP 


ABSORBED 


REJECTED 


POWER INPUT 


COP 


(degrees F) 


(MBtu/hr ) 


(MBtu /hr ) 


(kw) 


(MBtu/hr ) 




40 


11.5 


18.0 


1.80 


6. 14 


2.93 


50 


13.3 


20. 7 


2.06 


7.03 


2.94 


60 


15.2 


23.2 


2.30 


7.85 


2.96 


70 


17. 1 


25.9 


2.53 


8.63 


3.00 


80 


19. 0 


28.4 


2.76 


9.42 


3.01 


90 


20.9 


31. 1 


2. 99 


10.20 


3.05 


100 


22.8 


33.7 


3. 22 


10.99 


3.07 


1 10 


24.7 


36.3 


3. 45 


11.77 


3. 08 



2.3 Operation of System 



The solar assisted heat pump system is designed to 
supply the complete heating requirements of the residence 
with minimal energy consumption. The required heat is 
provided by one of three modes. These modes are direct 
heating from storage, heat pump heating and auxiliary 
electric resistance heating. 



The collector loop uses a 40Z ethylene glycol-water 
mixture to prevent freezing during low ambient temperature 
conditions* The mixture is pumped through the collector/ 
heat exchanger loop at a rate of .025 gpm/sqft of collector. 
The pump Is controlled by a control mechanism which turns 
the pump off when the temperature of the fluid leaving the 
collector Is less than the storage tank temperature. The 
pump resumes operation when the temperature of the fluid In 
the collector reaches the storage tank temperature plus 10 
degrees F. This control sequence prevents the loss of 
thermal energy from the tank by reradlatlon through the 
collector at night or cloudy days when solar Insolation 
decreases to a minimum. 

When the house thermostat calls for heat and the 
storage tank temperature Is greater than or equal to 110 
degrees F water from the tank Is pumped through a heating 
coll in the house air duct while the house air fan 
circulates warm air through the house. This is direct 
heating from the storage tank. If heating is required and 
the storage tank temperature Is below 110 degrees F the 
direct heating mode is disabled and the control valves 
circulate the water through the heat pump heat exchanger. 
The heat pump now provides heat to the house. On occasion 
the capacity of the heat pump will be unable to meet the 
house heating requirement In which case the built-in house 
air stream electric resistance heater of the heat pump will 



provide the additional required heat* If the storage tank 
temperature falls below 40 degrees F and heating is still 
required then the heat pump is shut down and all the 
required heating is provided by the electric resistance 
h ea t e r • 

If the ambient air temperature exceeds 65 degrees F the 
heating requirements decrease to zero and all heating modes 
become inoperative. 

2.4 Computer Model 

The solar assisted heat pump system shown in Figure 2 
was modeled using the Transient Simulation Program (TRNSYS) 
developed at the University of Wisconsin. The program 
allows one to make a transient (time dependent) simulation 
of the overall system by interconnection of individual 
component models. The TRNSYS modeling components used are 
shown in Figure 4 which is a flow diagram of the computer 
model . 



The reader component is used to input hourly weather 
data. The simulation is run for one month of the heating 
season. Only ambient temperature and solar insolation 
(measured on a horizontal surface) data are needed for the 
simulation. Actual weather data obtained from the 
Pennsylvania State University weather station for the months 
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FIGURE 4 — TRNSYS computer model flow diagram (Terms 

are defined in Appendix B) 
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of December 1960 and December 1972 are used. For the past 
thirty years of available weather data, these years 
represent the highest and the lowest values for December 
solar Insolation, respectively, for State College, Pa. 

The radiation processor component converts the solar 
Insolation or radiation on a horizontal surface (H), to 
solar radiation on a tilted surface (UT). The collector Is 
tilted at an angle of 50.5 degrees. By using Mode 3 of the 
radiation processor component the total radiation on a 
tilted surface Is a combination of the beam radiation, 
diffuse radiation and reflected radiation. The radiation 
processor calculates these values from the horizontal 
surface radiation, the collector orientation and the ground 
reflectance which Is .2 for bare ground. 

The liquid collector-storage module Is a single 
component used to model the solar collector, ethylene 
glycol/water heat exchanger, storage tank and associated 
pumps and relief valve. A schematic of the module Is shown 
In Figure 5. The combining of these components Into one 
subsystem reduces computation time and expense. The 
collector fluid, 40% ethylene glycol solution. Is pumped at 
the collector manufacturers recommended rate of .025 
gpm/sqft of collector. Across the heat exchanger the 
capacity rate ratio Is held at 1.0 [l.e. (mCp)eg (n;)Cp)w ]. 



m 




The initial storage tank temperature at the beginning of 
December is assumed to be 70 degrees F. It turns out that 
any initial temperature in this range gives nearly identical 
results over the one month period. The relief valve is set 
to maintain the tank temperature below 212 degrees F. 




FIGURE 5 - TRNSYS liquid co 1 lec tor /s tor age module 



The water-to-alr heat pump component uses actual heat 
pump performance characteristics. The model is q ua s ;l -s t ea d y 
state in nature. That is, at any instant in the "transient** 
simulation, the heat pump performs as. if it were at steady 
state, having for the moment constant source and sink 
temperatures and hence a unique heat addition, heat 
rejection and work input. The heat pump is modeled to 
operate in the heating mode only. A parallel auxiliary 



electric heating element is included to add Just enough 



energy to meet the heating load when the heat pump capacity 
is less than the load or when the storage tank temperature 
falls below 40 degrees F* Figure 6 is a schematic of the 
water-to-air heat pump component. 



ROOM 

AIR 




FIGURE 6 - TRNSYS water-to-air heat pump component 



The heating load component determines the heating load 

of the residence by use of the energy /degree day concept. 

The concept is extended to estimate the hour-by-hour heating 

load by multiplying the heating requirement of the residence 

< 

(energy /degree-hour ) by the difference between the required 
room temperature and the outdoor ambient temperature. 



The quantity Integrator simply keeps a continuous total 
of various system outputs throughout the simulation* It 
works slmlllar to a KWH meter which continuously to<:als the 
amount of electrical energy consumed. 

The printer components control the display of the 
system outputs. In this simulation some of the outputs are 
printed at hourly Intervals and some at dally Intervals. 

Each of the TRNSYS components requires a set of 
defining parameters and Inputs to produce an output. The 
parameters describe the operating characteristics of the 
system component being modeled. Many of these parameters 
have been discussed previously. A complete list and 
description of the typical parameters used In this 
simulation for each component can be found In Appendix B. 
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CHAPTER 3 



3«1 Presentation and Discussion of Results 

The collector efficiency, which is a measure of the 
useful solar energy being transferred to the heating system, 
is shown in Figures 7 and 8 for various values of collector 
area and storage capacity. The collector efficiency is 
maximum where the collector area is small and the storage 
capacity is large. This can be shown also by the collector 
efficiency equation when one realizes that the inlet fluid 
temperature to the collector (from storage) will be a 
minimum when the ratio of collector area to storage capacity 
is a minimum. In other words looking at equation (1) again, 

- .69-.7I (Tfl-Ta)/HT] , (1) 

low values of Tfl (resulting from a small collector area to 
storage capacity ratio) will reduce the second term on the 
right hand side resulting in a larger value for the 
collector efficiency. The collector efficiency increases 
for increases in the storage capacity and decreases in the 
collector area since both of these trends decrease the 
temperature of the fluid entering the collector. Figure 7 
also shows that increasing the storage capacity beyond about 
2 gal/sqft of collector area is of little benefit since the 
collector efficiency remains relatively constant. 
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FIGURE 7 — Collector efficiency vs storage capacity 
to collector area for solar assisted heat 
pump system in State College, Pa. 
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FIGURE 8 - Collector efficiency vs collector area for 
solar assisted heat pump system in 
State College, Pa. 
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The effect of weather on the collector efficiency can 
be seen by comparing the curves for 1960 and 1972 in Figures 
7 and 8. The 1972, or minimum solar insolation, curves 

reflect higher collector efficiencies for identical system 
parameters (collector area and storage capacity). This 
result is not as easy to predict from equation (1) since all 
three variables on the right side change. From the 
graphical results, however, it can be concluded that for 
poor weather the predominant effect on the efficiency 
equation is a decrease of the second term on the right hand 
side. This implies that the magnitude of the drop in 
entering fluid temperature must be larger than the drop in 
solar insolation, indicating relatively low entering fluid 
temperatures (storage tank temperatures) would be expected 
for poor weather conditions. Figure 7 also shows that 
weather extremes have a negligible effect on the collector 
efficiency of systems with small collector area and a 
significant effect on the collector efficiency of systems 
with large collector area. 

Values of the System Coefficient of Performance (SCOP) 
for various collector area and storage capacity 

configurations are shown in Figures 9 and 10. The SCOP is 
defined by the following equation: 

total heating load 

SCOP (7) 

electric energy used by heat pump & aux heater 
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SYSTEM COEFFICIENT OF PERFORHANCE 




FIGURE 9 - System coefficient of performance vs storage 
capacity to collector area for solar assisted 
heat pump system In State College, Pa. 
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figure 10 - System coefficient of performance vs collector 

area for solar assisted heat pump system In 
State College, Pa. 
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Following conventional practice, the energy required for 
fans and pumps is considered negligible in comparison with 
other energy inputs* Since the total heating load 
requirement for the residence is independent of the heating 
system configuration, the SCOP will depend only upon the 
amount of energy used to run the heat pump and the energy 
required for the auxiliary heater* The required energy 
input from these two sources is reduced, however, when there 
is a contribution from direct heating. 

In Figures 9 and 10 it can be seen that all the curves 
approach a SCOP of 3.0 as the overall system size (collector 
area and storage capacity) increases. This is because as 
the storage capacity to collector area ratio Increases the 
system becomes more and more dependent upon the heat pump 
alone (l.e. the direct heating contribution and the 
auxiliary heat contribution to the heating load decrease). 
Specifically, the curves converge to a value of 3.0 since 
that is the average heating COP of the AAF heat pump over 
its range of evaporator inlet (storage tank) temperatures* 

Figures 9 and 10 both show that the maximum SCOP is 
obtained with a large collector area and a small storage 
capacity. The factor responsible for driving the SCOP past 
3.0 (the average heating COP of the heat pump) is the direct 
heating contribution. Without a direct heating contribution 
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the SCOP will be llinlted to a maximum value equal to the 
heating COP of the heat pump. Also from Figure 10 it is 
seen that for a SCOP greater than 3.0 the solar assisted 
heat pump system must have as a minimum a collector area 
greater than about 650 sqft. 

The collector area of 650 sqft in Figure 10 is also 
important in that it is the approximate collector area for 
which all the different storage capacity curves resulting 
from 1960 weather data intersect each other at a SCOP equal 
to 3.0. The significance of this intersection is that it 
represents the point at which the total residence heat 
requirement is provided by the heat pump and system 
performance is independent of the storage capacity. This 
means that at 650 sqft of collector area the energy input 
from the solar collector array is equal to the energy 
withdrawn from the storage tank by the heat pump, resulting 
in neither the benefit of a direct heating contribution or 
the requirement for auxiliary heat. By ex t r ap ola t ioi^ it is 
estimated that a similiar Intersection point will occur for 
the storage capacity curves resulting from 1972 weather data 
at about 1800 sqft. Since 1960 and 1972 represent the 
extremes of solar insolation it is logical to expect that 
for the average year the intersection point will be 
somewhere between 650 and 1800 sqft, say 1200 sqft. Thus, 
based on a year of average weather the solar assisted heat 





4 







pump system for State College would require a collector area 
of 1200 sqft or greater to obtain a SCOP greater than 3.0 
(the heat pump heating COP). 

For collector areas below the intersection point the 
residence heat load is supplied essentially by the heat pump 
and auxiliary heat and the SCOP increases as the storage 
capacity increases. For collector areas above the 
intersection point the residence heat load is supplied by 
the heat pump, auxiliary heat and direct heat and the SCOP 
decreases as the storage capacity increases. 

This situation suggests that a two tank storage system 
would be beneficial. A smaller high temperature tank would 
be charged first by the solar collectors and be used for 
direct heating. A second larger tank in the medium 
temperature range would be charged with the excess solar 
energy and be used as a heat source for the heat pump. This 
combination would allow the direct heating contribution to 
be maximized. Whether or not the advantages of a two tank 
system outweigh the added cost and complexity would have to 
be studied in more detail. Figures 9 and 10 also indicate 
that for a single storage tank system the SCOP is more 
strongly influenced by the collector area than by the 
storage capacity, since in general the SCOP increases more 
for an increase in collector area than for an increase in 
storage capacity. 
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Comparing the curves resulting from 1960 weather data 
(good) and 1972 weather data (poor) in Figures 9 and 10 
shows that the weather data used for a solar assisted heat 
pump system simulation has a major effect on the resulting 
system performance characteristics. The SCOP obtained for a 
given collector area-storage capacity configuration can vary 
by more than 250Z. If non-representative weather data are 
used in the design calculations, a system design could 
result which is inadequate for one situation or overdesigned 
for another. 

Figures 11, 12 and 13 show the precentage contribution 
made by the heat pump, direct heat and auxiliary heat to the 
total heating load for a range of collector areas and 
storage capacities. For a given collector area the heat 
pump contribution Increases as the storage capacity 
increases. Accordingly, the percentage of the total heat 
load supplied by auxiliary heat decreases. Looking at the 
bar graphs for 1960 and 1972 (representing weather extremes) 
it is seen that poor weather, as intuitively expected, 
increases the auxiliary heat requirements. The amount of 
increase is most severe around collector areas of 600 sqft 
where it averages about 403!. The percentage contribution 
made by the heat pump, for a fixed storage capacity to 
collector area ratio, peaks in the middle range of collector 
areas. As predicted previously by Figures 9 and 10, the 
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FIGURE 11 - Source of heat for solar assisted heat 

pump system in State College, Pa* 
at 200 sqft of collector area 
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FIGURE 12 — Source of heat for solar assisted heat 
pump system In State College, Pa. 
at 600 sqft of collector area 



- 31 - 













• :• ••••*• 




HEAT POMP 




AUX BEAT 





DIR^ HKAT 



COLLECTOR AREA - 1000 SQPT 




STORAGE CAPACITY TO COLLECTOR AREA 



FIGURE 13 - Source of heat for solar assisted heat 
pump system in State College, Pa. 
at 1000 sqft of collector area 
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heat pump contribution is near 100% for a collector area 
near 650 sqft (for 1960 weather data)* At large collector 
areas (1000 sqft. Figure 13) the heat pump contribution is 
smaller for good weather than for poor weather because the 
direct heat contribution lessens the heat requirement to be 
provided by the heat pump* 

Direct heating becomes a significant contribution when 
the collector area is large and the storage capacity is 
small because the large amount of collected solar energy 
only has to heat a small amount of water, resulting in a 
high temperature level in the storage tank* Decreasing the 
collector area causes the percentage contribution of direct 
heating to decrease rapidly so that for small collector 
areas the direct heating contribution is non-existent* For 
a sufficiently large collector area, increasing the storage 
capacity will decrease the direct heating contribution* 
Also, poor weather will decrease the likelihood of any 
direct heating contribution* In Figure 13, for a storage 
capacity of 1 gal/sqft of collector, the direct heat 
contribution for poor weather (1972) is only 5% compared to 
A6Z for good weather (I960)* As indicated by previous 
figures, the bar graphs of Figures 11, 12 and 13 show that 
the best performimg configuration of the solar assisted heat 
pump system is a large collector area with a small storage 
capacity, mainly because of the advantage of more direct 
h ea t ing * 
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Table 2 lists the energy usage of a system with 1000 
sqft of collector area and a storage capacity of 1 gal/sqft 
of collector* Energy data is given for both "good** and 
**poor** weather* It is clear from the percent deviations 
listed that weather has a major impact on the system 
performance and energy requirements* Note that when the 
direct heating contribution falls off during poor weather 
the majority of the difference is made up by auxiliary heat* 
Several comparisons can be made from Table 2 besides the 
deviations shown* For instance, assuming an electric rate 
of $*035/kwhr for State College and the system configuration 
given in Table 2, it would cost between $28 and $73 
(depending upon the weather) to heat the house for the month 
of December* 



TABLE 2 - A one month comparison of system performance 
characteristics for a solar assisted heat pump 
system in State College, Pa* ; Collector area - 
1000 sqft and water storage capacity - 1000 gal 



Collector Efficiency (%) 

Sys Coefficient of Performance 
Total Energy Collected (kwhr) 

Total Heating Load (kwhr) 

Heat Pump Contribution (kwhr) 
Direct Heat Cont r ibu t ion (kwh r ) 
Aux Heat Contribution (kwhr) 

Energy Input to Heat Pump (kwhr) 
Energy Input to Aux Heat (kwhr) 



c 1960 


Dec 1972 


Percent 


athe r 


Wea the r 


Dev ia t ion 


Data 


Data 


From 1960 


36 


44 


22.2 


5. 14 


1.90 


-6 3.0 


3648 


1778 


-5 1.2 


4161 


3950 


-5.1 


2151 


2508 


16.6 


1902 


193 


-89.9 


100 


1242 


1142 


709 


84 1 


18.6 


100 


1242 


1 142 
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Since the collector array Is mounted on the roof of the 
residence, in reality an upper limit of collector array 
size is set by the area of the roof* Taking this into 
consideration a collector area of 600 sqft was used for the 
parametric studies of room thermostat setting and ethylene 
glycol/water heat exchanger size* Figure 14 shows for a 
collector area of 600 sqft the effect on SCOP of lowering 
the room temperature to 65 degrees F as the storage capacity 
is varied* It is interesting to note that the SCOP is 
increased when the room temperature is maintained at 65 
degrees F versus 70 degrees F* Recalling that the SCOP is 
defined by equation (7) as 

total heating load 

SCOP (7) 

electric energy used by heat pump & aux heater 

one might expect that since the heating load is reduced by 
about 13% the SCOP would decrease, but the reduction in 
energy input to the system is even greater* The heat pump 
input energy is reduced by 13% and the auxiliary heat input 
is reduced to nearly zero resulting in a net Increase in the 
SCOP* As shown by the nearly horizontal curves of Figure 
14, and as previously stated for a collector area of 600 
sqft, the storage capacity size has little effect on the 
SCOP* 



Table 3 shows how lowering the room temperature to 65 
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STORAGE CAPACITY TO COLLECTOR AREA (GAl/SQFT) 



FIGURE 14 - Effect of room temperature on system coefficient 

of performance for a solar assisted heat 
pump system in State College, Pa« 
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TABLE 3 - Comparison of energy requirements for room 
temperatures of 65 F and 70 F for a solar assisted 
heat pump system In State College, Pa. 



DEC 1960 WEATHER DATA 


Collector Area 


“ 600 sqft 




Storage Capacity 


« 6000 gal 




Room Temp 


Room Temp 




65 F 


70 F 


Total Heating Load (kwhr) 


360A 


4161 


Heat Pump Contribution (kwhr) 


3598 


4120 


Direct Heat Contribution (kwhr) 


0 


0 


Aux Heat Contribution (kwhr) 


. 03 


34 


Energy Input to Heat Pump (kwhr) 


1193 


1374 


Energy Input to Aux Heat (kwhr) 


. 03 


34 



DEC 1972 WEATHER DATA 



Total Heating Load (kwhr) 

Heat Pump Contribution (kwhr) 
Direct Heat Contribution (kwhr) 
Aux Heat Contribution (kwhr) 

Energy Input to Heat Pump (kwhr) 
Energy Input to Aux Heat (kwhr) 



Collector Area - 600 sqft 
Storage Capacity ■ 6000 gal 



Room 


Temp 


Room 


Temp 


65 


F 


70 


F 



3393 3950 



2622 


2683 


0 


0 


765 


1263 


888 


909 


765 


1263 



degrees F affects the energy inputs of the system. A 
storage capacity of 6000 gallons was chosen since Figure 10 
shows that for 600 sqft of collector area the maximum SCOP 
is obtained with a storage capacity of 10 gal/sqft of 
collector area. The economics of lowering the room 
thermostat setting can be obtained from the values in Table 
3. Again assuming an electric rate of $.035/kwhr for State 
College it is seen that for the system configuration given 
in Table 3 approximately $8 to $18 can be saved (depending 
upon the weather) for the month of December by lowering the 
thermostat setting 5 degrees. 

Tables 4 and 5 are the system performance 
characteristics for good (1960) and poor (1972) weather 
resulting from variations in the ethylene glycol/water heat 
exchanger effectiveness. If the capacity rate ratio is held 
constant at a value of 1.0 (i.e. Cmin/Cmax - 1.0) the heat 
exchanger effectiveness (E) is defined by the following 
equa t ion : 



E - NTU/(1+NTU) (8) 

where NTU * Number of transfer units 

(size parameter) 

Thus by holding the capacity rate ratio constant at a value 
of 1.0, variations in the heat exchanger effectiveness (E) 
actually represent variations in the size of the heat 



-3 8- 



TABLE A 



- System performance factors of solar assisted heat 
pump system for varying values of heat exchanger 
effectiveness based on DEC 1960 weather data for 
State College, Pa. 



STORAGE COLLECTOR CONTRIBUTION TO HEAT LOAD (Z) 



HXGR 

EFF 


TANK 
TEMP (F) 


EFFICIENCY SYSTEM 
(Z) COP 


HEAT 

PUMP 


DIRECT 

HEAT 


AUX 

HEAT 


(Collector 


Area-600s qf t Storage Capacl ty »600gal ) 


. 6 


115.8 


4A . 3 


2.68 


73 


lA 


13 


. 7 


117.3 


AA . 7 


2.71 


72 


15 


13 


.8 


118.3 


A5. 0 


2. 7A 


71 


16 


13 


(Collector 


Area-600sqf t Storage Capaclty*2 AOOga 1) 


. 6 


87. 1 


A7. 7 


2.93 


99 


0 


1 


. 7 


88.8 


A8.0 


2.9A 


99 


0 


1 


.8 


90.0 


A8. 3 


2.95 


99 


0 


1 


(Collector 


Area“600sqft Storage Cap ac Ity “6000ga 1 ) 


.6 


76.7 


A8. 5 


2.95 


99 


0 


1 


. 7 


78.0 


A9.0 


2.96 


99 


0 


1 


.8 


78.9 


A9. A 


2.96 


99 


0 


1 



TABLE 5 - System performance factors of solar assisted heat 
pump system for varying values of heat exchanger 
effectiveness based on DEC 1972 weather data for 
State College, Pa. 



STORAGE COLLECTOR CONTRIBUTION TO HEAT LOAD (Z) 



HXGR 

EFF 


TANK 
TEMP (F) 


EFFICIENCY SYSTEM 
(Z) COP 


HEAT 

PUMP 


DIRECT 

HEAT 


AUX 

HEAT 


(Collector 


Area“600s qf t Storage Capacity“600gal ) 


. 6 


AO. 2 


A8.8 


1 . A7 


A6 


1 


53 


. 7 


39.8 


A9. A 


1.A8 


A6 


1 


53 


.8 


AO . 1 


A9. 8 


1 . A8 


A 7 


1 


52 


(Colle ctor 


Area*600sqft Storage Capacity-2A00ga 1) 


.6 


AO.O 


52.6 


1 .62 


57 


0 


A3 


.7 


AO.O 


53.3 


1.63 


58 


0 


A2 


.8 


AO.O 


53.8 


1 .63 


58 


0 


A2 


(Co lie ct or 


Area-600sqft Storage Capacity“6000ga 1) 


. 6 


AO.O 


52.0 


1.81 


67 


0 


33 


. 7 


AO.O 


52.7 


1.82 


68 


0 


32 


.8 


AO.O 


53. 3 


1.83 


68 


0 


32 
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exchanger* The smaller the value of £ the smaller Is the 
size of the heat exchanger* As can be seen from Tables 4 
and 5 varying the heat exchanger effectiveness (l*e* size) 
of the ethylene glycol/water heat exchanger has negligible 
effect on the performance of the system* The SCOP Increases 
very slightly as the heat exchanger size Increases 
regardless of the weather for the effectiveness range 
studied* 
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CHAPTER 4 



4 • 1 Cone lus ions 

The maximum collector efficiencies are obtained for 
system configurations with small collector areas and large 
storage capacities because the storage tank temperature and 
thus the temperature of the fluid entering the collector is 
kept low. Similarly, the collector efficiency of a given 
system is improved during periods of poor solar insolation 
due to the reduced entering fluid temperature. The 
collector efficiency increases as the storage capacity 
increases and as the collector area decreases since the 
combination of these result in a lower temperature of the 
fluid entering the collector. Weather extremes have little 
affect on the collector efficiency of systems with small 
collector area, however, the affect on the collector 
efficiency is significant for systems with large collector 
area. The study results also show that for the solar 
assisted heat pump system a storage capacity beyond 2 
gal/sqft of collector area is of little benefit since the 
collector efficiency remains relatively constant. 

For a solar assisted heat pump system the SCOP is 
limited to the heating COP of the heat pump unless the 
combination of good solar insolation and adequate collector 



area result in a heating contribution from direct heating. 
The SCOP will be largest for systems with large collector 
areas and small storage capacities. There Is a unique 
collector area at which the SCOP of the solar assisted heat 
pump system will equal the heating COP of the heat pump 
(l.e. lOOZ of the total residence heating requirement Is 
provided by the heat pump) for a wide range of storage 
capacities. This unique collector area Is dependent upon 
the solar Insolation received and will vary In value 
according to the weather. For the solar assisted heat pump 
system modeled for State College In this study It Is 
estimated that about 1200 sqft of collector area would be 
required to provide lOOZ heat pump heating under average 
December weather conditions. 

The SCOP of the solar assisted heat pump system 
Increases with Increasing collector area and In general Is 
more strongly Influenced by the collector area than by the 
storage capacity. For a system with a collector area 
greater than that resulting In 100% heat pump heating, the 
residence heating load Is supplied by contributions from the 
heat pump, auxiliary heat and direct heat and the SCOP 
decreases as the storage capacity Increases. For a system 
with a collector area smaller than that resulting in 100% 
heat pump heating, the residence heating load Is supplied 
essentially by the heat pump and auxilary heater and the 



SCOP increases as the storage capacity increases. This 
situation suggests the use of a two tank storage system with 
a high temperature tank being used to maximize the direct 
heat contribution and a medium temperature tank for use as 
the heat pump heat source* 

For a given collector area the percentage contribution 
of the heat pump increases as the storage capacity increases 
whereas the percentage contribution of the auxiliary heater 
decreases* Direct heating becomes a significant 
contribution to the total heating load for solar assisted 
heat pump systems with large collector areas and small 
storage capacities since the storage is heated to a high 
temperature level* Decreasing the collector area causes the 
percentage contribution of direct heating to decrease 
rapidly so that for small collector areas the direct heating 
contribution is non-existent* 

When designing or simulating solar assisted heat pump 
systems it is essential that accurate and realistic weather 
data be used* The weather data used has a significant 
effect on the performance characteristics of the system* 
For the system modeled in this study the SCOP varied by as 
much as 250% for some co 1 le ct o r /s t o r ag e configurations when 
"poor” versus "good" weather data were used* The percentage 
contributions of the heat pump, auxiliary heater and direct 
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heating to the total heating load are similarly affected as 
was shown In Figure 13. A simple economic analysis of a 
system In State College with 1000 sqft of collector and 1000 
gallons of storage shows that monthly heating costs can vary 
from $28 to $73 depending upon the weather. 

Federal energy conservation guidelines call for 
reducing thermostat settings to 65 degrees F for winter 
heating. This study found that reducing the room 
temperature from 70 degrees F to 65 degrees F Increased the 
SCOP approximately 5% and reduced significantly the amount 
of auxiliary heat required. For a system configuration of 
600 sqft of collector area and 6000 gallons of storage 
capacity a simple cost analysis shows that In State College 
approximately $8 to $18 per month (depending upon the 
weather) can be saved by lowering the thermostat setting 5 
degrees . 

The performance of the solar assisted heat pump system 
Is not significantly affected by the effectiveness (l.e. 
size) of the ethylene glycol/ water heat exchanger. Since 
this Is the case. It appears that the heat exchanger could 
be of relatively simple design and low cost. 

When reviewing the results presented In this study It 
must be remembered that several assumptions have been made 



in modeling and simulating the operation of the residential 
solar assisted heat pump system. The majority of these 
assumptions are inherent in the TRNSYS components but others 
are externally controlled. For instance the effect of wind 
velocity has been assumed negligible which in a **real world” 
installation could reduce the solar energy collected by the 
collector. Also a ground reflectance value for bare ground 
has been assumed. In a "real world" installation the 
reflectance value might be higher if the ground was snow 
covered resulting in an increase of reflected solar 
radiation received by the collector. Overall, however, it 
is expected that the model performance will be somewhat more 
optimistic than that of a "real world" installation. 



Although domestic hot water heating was not included in 
the solar assisted heat pump model of the residential 
heating system, the collector array sized for the solar 
assisted heat pump system would most likely provide all of 
the domestic hot water requirements in the summer months. 
During the heating season the system would be able to 
provide preheating of the domestic hot water much of the 
t ime . 



At first glance the TRNSYS modeling program seems very 
complicated but actually the building block modeling 
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technique Is fairly straightforward. Individual system 
components (i.e. pumps, solar collectors, storage tank, 
etc.) are described by component models which are connected 
together to form a complete system. It Is better to keep 
the model as simple as possible since problems, such as 
failure of the differential equation solution techniques to 
converge within the component models, are likely to result 
if the model Is too elaborate. Some modes of the TRNSYS 
components use empirical equations which have been developed 
in SI units requiring the entire simulation to be run in SI 
units if these component modes are used. If these component 
modes can be avoided in the model the more familiar English 
units can be used. The TRNSYS simulations require hourly 
weather data as input. The assembly and handling of such 
data can become quite cumbersome if punched cards are being 
used, especially if the simulation is run for a full year. 
This problem can be alleviated by using a tape input. The 
computer time required for running a full year simulation of 
a solar assisted heat pump system (based on this one month 
run) will be about 200 seconds and cost approximately $60 if 
output is printed for each hour. 

4.2 Recommendations 

To obtain a complete picture of the solar assisted heat 
pump system, a detailed economic analysis of the system 



- 46 - 



alternatives should be performed. 



This would determine the 



most cost effective combination of such things as collector 
area and storage capacity. An economic analysis would also 
be an effective method to compare the solar assisted heat 
pump system with conventional heating systems. 

An entire year or at least a complete heating season of 
weather data should be used when simulating the solar 
assisted heat pump system rather than a single month. This 
will give a more realistic picture of the system's seasonal 
performance characteristics. Limiting the weather data to 
the heating season will significantly decrease the amount of 
data and reduce computer costs. It is important however 
that the weather data be carefully selected to provide a 
realistic representation of the location being studied. 

Future studies on the solar assisted heat pump system 
should explore the feasibility and operating characteristics 
of a two tank storage system for optimizing the output of 
the solar assisted heat pump system. Also the results 
presented in this report and any future simulations using 
the TRNSYS program should be compared with other modeling 
techniques such as those being developed by the 
Architectural Engineering Department of the Pennsylvania 
State University. 
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BUILDING CRITERIA 



Wall: 



Glass: 



Doors: 



Ceiling: 



Indoor Design Temperature ■ 70F 
Outdoor Design Temperature » 7*F 



CONSTRUCTION DATA 



4" Face Brick 

3/4" Air Gap (Mean Temp. = 50, Temp. Diff. ■ 10, £ ■ .82) 
Building Membrane (Permeable Felt) 

1/2" Plywood (Douglas Fir) 

3 1/2" Fiberglass Insulation 
1/2" Gypsum Wall Board 

2x4 Stud Wall Construction, 16" O. C. 



Area is given on plan; e.g. Ac 15 
All Double Pane 

(1/4" Air Gap, Wood Sash, 80% Glass) 
3' X 7' 

1 1/2" Solid Wood 

Metal Storm Door (Front Door only) 

6" Fiberglass Insulation 
1/2" Gypsum Board 

2x6 Ceiling Joists, 16" O. C. 



I 
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U-Valuea; 



Walls: 



R 



1977 Handbook of 
Fundamentals Reference 
Page No. Table No. 



Outside Air Film (15 mph wind) 


./7 


7.2. \\ 


4" Face Brick 




^ZZiLS^ 



3/4" Air Gap /.o7 ZX.IZ- 

(Mean Temp, = 50, Temp. Dlff.® 10, 

E = . 82) 



Building Membrane (Permeable Felt) 


• i 


22, 


1/2" Plywood (Douglas Fir) 


r 

1 i 


22.73 


3 1/2" Fiberglass Insulation 


UJ.OO 


Z2.14 


1/2" Gypsum Wall Board 






Inside Air Film 




ZZ.lf 


TOTAL RESISTANCE 


I4.i3 




FINAL U-VALUE 






Class : 


u 




Double Pane 
(1/4" Air Gap) 




22,74- 


Adjustment for Wood Sash, 80% Glass 
(Ux adjustment) 




22, z4 


FINAL'U-VALUE 


.55' 




Door; 


u 




1 1/2" Solid Wood 
with metal storm door 


.33 


IZ.2S 


FINAL U-VALUE 


.33 




Ceiling : 


R 




Air Film (Still Air, E = . 90) 


,6./ 


7ZAL_ 


6" Fiberglass Insulation 


ZZ.oo 


zz,i4 


1/2" Gypsum Wall Board 


.V5- 


22 JZ 


Inside Air Film 


,6/ 


^^,11 


TOTAL RESISTANCE 


zi.^n 




FINAL U-VALUE 


.c4z. 




Kitchen Wall i^^ainst Garaj^e: 


R 




Air Film (Still Air, E , 90) 




ZZJl 


3 1/2" Fiberglass Insulation 


11,00 


zz,i4 


1/2" Gypsum Wall Board 




ZZJ2 


Inside Air Film 




zz,n 


TOTAL RESISTANCE 


/Z«/ 




FINAL U-VALUE 


,oi% 





/-/) 

ZlA 

2L 













!-f\ 



IvT 



7^ 
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HEAT LOSS SUMMARY 



Location (City ) 

Design IDBT » 70F 
Design ODBT = 7^ F 



Design Temp. Diif. « 



63 



Type of BuUdin g 2 STZWy k€Sit&4Tlf\L 



STRUCTURAL 

COMPONENT 

Rloo^ 


AREA. 

^SG^FT) 


U-VALUE 


UA 


heat LOSS, 
Btuh 


lV/)ZiS 


m 


.0(ol 


G^.oo 




GifiSS, 


141 


.5X 


77-S5' 


4S2(a 


FLoo^ 


A ssuAfe- 


a/O A/eAT los‘^ 




I)oo/^S 


Zl 
2/ . 


/Azyj 


/OSS' 






.on 


%.4Z 


S^o 


►A/b R.^/^ 


UJ/lZ-Z-S 


im 


.okS 


H-Az 


4^11 . 


(TiT/Jss . 


ill 


.5^ 


6l.cs- 


3S4& 


C£ILIa/C 


iios 


.OAZ 


/I6.41 


ziz4 






SUB 


TOTAL 


ZZZ43 



/fSSUM, 



% /-OO 



/7S^7 rr^ 



INFILTRATION: No. of Air Changes /44Cfr/l \ll/t4K 

Q^UQCF/A)AT 

POSITIVE VENTILATION: 0 CFM 



7^77? 



TOTAL 



S2Z2Z. 






HEAT LOSS 
7'OTO'DbT" 



32222. 



Btuh ^ — y y 

^F ind-out temp dilf ““ ^ • f 



10- 7 
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COOLING LOAD CALCULATION 



SOURCES; Chap. 22, 1972 ASHRAE Handbook of Fundamentals 
Chap. 25 , 1977 ASHRAE Handbook of Fundamgntals . 

Location (City) St^te’ CoLcece', P/q 
Design Conditions 75 I DBT, S~7 O DBT 

Daily Range ^3 F 

Design Equivalent Temperature Differences: 

Brick Walls g.3 F CeiUng Roofs F 

Wood Doors /S^ ^ F Floors O F 



Design Window Factors: 

North /3 Btuh/ft^ 
-<^3 Btuh/ft^ 



Wes t Btuh/ft^ 

Sout h Btuh/ft^ 



COOLING LOAD SUMIvlARY 



Walls 11^^ 

Windows /(o 

Doors ?jZ.I 

Ceiling ! (ol I 

Floor O 

Occupants 

Cooking . 12.00 

ML 7 Infiltration 2/ /4 

Total /3>Zo4 



Btuh 



Btuh, Total Sensible Heat (T£ 



Total Heat ■ 1.3 (TSH) « Btuh 



Equipment Size 



Total Heat, Btuh 
1 2, 000 Btuh/ton 



. /. 44 - 
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tons 



APPENDIX B 
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DEFINITION OF TERMS USED IN TRNSYS COMPUTER MODEL FLOW DIAGRAM 



DELTAU 


- Total change in storage tank internal energy 


H 


- Solar radiation (total) on a horizontal surface 


Hb 


- Solar radiation (beam) on a horizontal surface 


Hd 


- Solar radiation (diffuse) on a horizontal surface 


HT 


- Solar radiation (total) on a tilted surface 


HB 


- Solar radiation (beam) on a tilted surface 


HD 


- Solar radiation (diffuse) on a tilted surface 


HI 


- Fraction of tlmestep that heat pump operated in 
water-air heating mode 


H2 


- Fraction of timestep that heat pump operated in 
air-air heating mode 


H3 


- Fraction of timestep that heat pump operated in 
cooling mode 


IGAM 


- Heat pump master control Integer; 0*off, 1-on 


Ml 


- Flowrate of water entering heat pump evaporator 


Mo 


- Flowrate of water returning from heat pump 
evaporator 


MODE 


- Mode of heat pump operation 


OPHRS - 


> Total time heat pump operated in heating mode 


QA 1 


‘ Instantaneous energy absorbed from liquid source in 
water-air heating mode 


QABS 


• Total energy absorbed by heat pump 


QAC 


• Instantaneous energy removed from load in cooling 
mode 


QAUXH - 


Auxiliary energy delivered to room in heating mode 


QDH 


Energy delivered to room in direct heating mode 


QDUMP - 


Rate at which energy is discarded by relief valve 


QENV 


Rate of energy loss from storage tank to 
surroundings 
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QL 


- Instantaneous room heating load 


QLOAD 


- Hourly heating load 


QRl 


- Instantaneous energy delivered to load in vater-air 
heating mode 


QR2 


- Instantaneous energy delivered to load in 4^ir-air 
heating mode 


QRC 


- Energy rejected to sink in cooling mode 


QREJ 


- Total energy delivered by heat pump 


QSHRTC 


- Cooling capacity required but not met by heat pump 


QTANK 


- Rate of energy delivery from storage tank to heat 
p ump 


Qu 


- Rate at which energy is transferred to storage tank 


QUSE 


- Hourly useful energy collected 


QW 


- Total work input to heat pump 


RAD 


- Hourly Incident solar radiation on collector 


T 


- Temperature of storage tank 


Ta 


- Ambient temperature 


TENV 


- Temperature of ground surrounding storage t,aink 


Ti 


• Temperature of water entering heat pump evaporator 


To 


• Temperature of water leaving heat pump evapprator 


TOTQL - 


• Total heating load 


TOTQU - 


• Total useful energy collected 


TOTRAD - 


■ Total incident solar radiation on collector 


TSTOR - 


■ Hourly storage tank temperature 


WAC 


> Instantaneous work input to heat pump in cooling 
mo d e 


WAH 


Instantaneous work input to heat pump in heating 
mode 


WORKIN - 


Total work input to heat pump in heating mode 


AE 


Change in internal energy of water in storage tank 



- 57 - 




"•'i 













I 




et - Angle of incidence of solar radiation on a tilted 

surface 







a f 















I 



TYPICAL TRNSYS COMPONENT PARAMETERS 



(A) UNIT 1, TYPE 9 READER 



Parameters Description Va lue 



1 N - number of values to be read 2 

2 At - time Interval at which data Is 1 hr 

pr ovl ded 

3 1 - position of value on card for 2 

which units are to be converted 

4 Ml - multiplication factor for that value 3*686 

5 A1 - addition factor for that value 0 



(parameters 3,4&5 convert Langleys to Btu/hr-sqft) 



6 L - logical unit number of data file 



(B) UNIT 2, TYPE 16 RADIATION PROCESSOR, MODE 3 



Pa rame ter 



Description 



Value 



1 

2 

3 

4 

5 

6 
7 




mode 

day of year at start of simulation 
latitude 

slope of collector surface 
orientation angle of collector 
sola r cons tant 
ground reflectance 



3 

335 

40.5 deg 

50.5 deg 
0 deg 

429 Btu/hr-sqft 

. 2 



(C) UNIT 3, TYPE 21 LIQUID COLLECTOR-STORAGE MODULE 



Parameter 



Descr Ip t Ion 



Value 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

1 1 
12 

13 

14 



A - collector area 

Me - collector fluid flowrate 

Cpc- collector fluid specific heat 

Ah - heat exchanger to tank flowrate 

Cph- specific heat of fluid In 
storage tank 

F' - collector efficiency factor 
~ t ransmlt tance-absorptance 
product 

UL - collector loss coefficient 
Tmax- temperature at which tank 
relief valve opens 
6 - c ol le c t or - 1 ank heat exchanger 

effect Ivenes s 

V - volume of storage tank 
P - tank fluid density 
U - tank energy loss coefficient 
r - ratio of tank height to tank 
dl ame ter 



600 8 qf t 
7500 Ib/hr 
.792 Btu/lb-F 
5940 Ib/hr 
1.0 Btu/lb-F 



.99 

.684 

.7 Btu /hr-s qf t-F 
212 deg F 

. 7 

80.1 cuf t 
62.4 Ib/cuft 
.05 Btu /hr-sqf t-F 
3 
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15 



To • Initial tank fluid temperature 70 deg F 



(D) UNIT 4, TYPE 20 WATER-TO-AIR HEAT PUMP 
Parameter Descr ip t ion 

1 CPI - specific heat of liquid 

entering heat exchanger A 

2 MDOTl - mass flow rate of liquid 

entering heat exchanger A 

3 TROOMH - constant room temperature 

for heating mode 

4 TROOMC • constant room temperature 

for cooling mode 

5 TMINl - minimum liquid source temp 

for water-to-air heat pump 

6 TMIN2 - minimum air source temp for 

air-to-air heat pump 

7 NDATAH - number of equally spaced 

heating data points 

8 NDATAC - number of equally spaced 

cooling data points 

9 ISAVEH - logical unit number of 

heating data file 

10 ISAVEC - logical unit number of 

cooling data file 

11 TCOOL - minimum ambient temp when 

cooling is allowed 

12 CAIR - MCp product for room air 

flow thru heat exchangers 
B and D 

13 EFF - effectiveness of heat 

exchanger B 

14 TSET - minimum liquid source temp 

for direct heating from 
liquid source 

15 ICOOL - cooling condenser selection 

integer, 0«heat exchanger A 

16 TREAT - maximum ambient temp when 

heating is allowed 



(E) UNIT 5, TYPE 12 HEATING LOAD, MODE 3 
Parameter Descr ip t ion 

1 mode 

2 UA - heating requirement of house 

(energy /deg re e hour) 

3 TR - room temperature 



Value 

1 . 0 Btu/lb-F 
2500 Ib/hr 
70 deg F 
77 deg F 
40 deg F 
1000 deg F 
8 
8 

31 

32 

1000 deg F 
648 Btu/hr-F 

.7 

1 10 deg F 
0 

65 deg F 



Value 

3 

-511 Btu/hr-F 
70 deg F 
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(F) UNIT 6, TYPE 24 QUANTITY INTEGRATOR 



No parameters required 



(C) UNIT 7, TYPE 25 PRINTER (HOURLY) 

Parameter De scr Ip t Ion 

1 Tp - time Interval at which printing 

Is to occur 



(H) UNIT 8, TYPE 25 PRINTER (24 hour) 



Parameter 



De scr Ip t Ion 



1 Tp - time Interval at which printing 

Is to occur 



Value 
1 hr 

Value 
24 hr 
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TRNSYS - A TRANSirNT SI'^UCATION P»fi';RAM 
rAOH TMC SOLAR fNEH&V lAB AT IHt JMVCRSIfY OF WlV(ON)lN. 
VERSION «• 1 V/? r/76 



SI’nula T ION 




0.0 




7.AA06 02 


2.5006-01 














TOL ERANCES 




- l.OOOE-02 


- l.OOOE-02 
















UNIT 


1 


TYPE 


U 


weather DATA READER 
















parameters 


6 






















2.000b 


)0 


1.C30E 


70 


2.030E 00 


3.666^ 00 


,0.3 


8.0006 


00 




- 




UNIT 


2 


TYPE 


1 6 


MAOIATION PROCESSOR, MODE 3 
















PAPAME TFRS 


7 






















3.000E 


30 


3. 350E 


02 


A.OSOE 01 


S.OSOfc 01 


0.0 


A. 2906 


02~ 


2.0706-01 






INPUTS 
If 2 


1 








» 












1 




0.0 
























UNIT 


3 


TYPE 


21 


liguiD COLLECTCR- 


STORAGE MOOULL 






' 










PAP AME TERS 


IS 






















6.3C0T 


02 


7. scot 


03 


7.9206-01 


5.9A06 03 


l.OOOE 00 


9.900E-01 


6. 0AOE-OI 


7 . 0006- 


01 




2. \20t 


02 


7.3006 


-01 


H.0lv7E 01 


0.2A3E 01 


^5.3 006-02 


3.0006 


00 


7.0 )06 01 








INPUTS 


5 
























2f 1 




If 1 




A, 1 


Af 2 


0, 0 














o 

o 




2.6C0E 


01 _ 


7.030E 01 


2.5006 03 _ 


5.0006 01 




- 








UNI T 


<• 


T YPP 


23 


mater 


TO AIR HEAT 


PUNP 
















PAP AMf T£0 S 


16 


















1 




1 .OJOE 


73 


2.500E 


03 


7.0306 01 


7.7006 01 


A. 0006 01 


' l.OOOE 


03 


8.0006 00 


4.030E 


00 




3. nOE 


01 


3.2036 


01 


1.0306 03 


O.A806 02 


7.7006-01 


1.1006 


02 


0.0 


4.5006 


01 




INPUTS 


. 5 
























3, 1 




3,2 




If 1 


5f 3 ^ 


■'0,0 ’ 














7.030E 


01 


2.5006 


03 


2. 8006 01 


-2.1A6E OA 


l.OOOE 00 


. 










UNIT 


S 


TYPE 


12 


HEATING LU\0<ENE RGY/OECR66 HOUR 


MODEL If MUOE 














PAR AHf TERS 


A 






















3.000E 


00 


-5.1 lOE 


02 / 


7.3306 01 


0.0 


p.o 


o 

0 

1 




o 

. 

o 


o 

. 

o 






INPUTS 


A 
























0, 0 




O 

o 




If 1 


Of 0 
















0.0 




0.0 




2. 8006 01 


0.0 














UNIT 


6 


TYPE 


2A Quantity integrator 
















INPUTS 


12 
























2f 1 




3f 3 




3f A 


if 5 


”3r> 


'Af 3 




Af A 


' ' Af 5 






Af 6 




A, 9 




A, 10 


5f 3 
















0.0 




0.0 




0.0 


0.0 


0.0 


o 

o 




o»o_ _ 


0*0 






0.0 




0.0 




0.0 


-2.iA6r^erA 















. UNIT 7 TYPC 25 PM I NTE R I hOURL Y I 

PAPAML'TERSl ■" 

l.OOOE 73 
INPUTS <» 

?• I yri’~ 3, 3 573 ^ 

RAO TSTOR UOSE OLOAO 

UNIT S TYRE ?5 P R I nTe R \ 2^^ mIJur] ' ’ 

parameters 1 

7.<»30E 31 

INPUTS 13 

3* ^ 6, I b, 2 6, 3 6f A 5 4f 6 4t 7 

_ 4,10 6|11 6 , 1 ? 

TOTRAO TUTUU UENV QTANK QOUMP WORR IN QAftS 




6 3 



I 




A 



T^>ANSfCNI SMULAI lUN 



SI A>» f I NO A I I I Mt • 0.0 
SfOPPINO AI TIME ■ y.AAOE 02 
T iMf SU P • 2.S00E-01 
OlfFf PENT I AL tuUAllCI feMPCH filUKANCF --I.OO0E-O2 
ALCtOKAIC CaNVEPGtNtC IfJLIPANtE •-1.000P-02 



RAO 



OOUMP 



I STOP 
y.OOOE 01 
WURK IN 

0.0 



oust 



0.0 



QABS 



0.0 



TIME 

ULOAO 

-2.l<»«»t 04^ 
ORE J 

0.0 



0.0 

OELU 


rniRAO 


TOTOU 


QENV 


qtana 


0.0 


0.0 


_ 0.0 


o 

• 

o 


0.0 


HRS 


OOH 


OAUX 


iotol 




0.0 


0.0 


1 

o 

• 

o 

1 

1 


-2.1AAE OA 





PAD 



0.0 



I IME 

TSfOrt 
6.750E 01 



l.OOOO 

OUSE 



0.0 



OLOAD 

-2.146E OS 



RAO 



0.0 



I IME 
TSTUR 
6 .A 66 E 01 



2.0000 

OUSE 



OLUAO 



0.0 



-2.146E OA 



PAO 



I 0.0 



I |Mfc 

fsruR 

6.179E 01 



3.0000 

UUSE 



0.0 



QLOAO 

-2.197E OA 



RAO 



0.0 



I IME 
TSfUR 
5.0OJE Cl 



A. 0000 
OUSE 



^ I 



0.0 



QtOAO 

-2.1A6E OA 



RAO 



0.0 



T IME 
f STOP 
S.609E 01 



5.0000 

OUSE 



0.0 



QLOAO 

-'2.197E OA 

























TIME - 


7A2.0000 












-X - 




RAO 


I STOP 


guSE 


giUAO 








* 






jO.O 


l.2/9f 02 


0.0 


1.3 HOE OA 


. 






- - 






1 


TIME • 


7Ai.D000 
















. PAO 


I SfUR 


QUSE 


QlUAD 














0.0 


1.2011 02 


0.0 


-1.329E OA 














1 


'l SfUK 




Tl ME 


• 7A<*.0300 












. PAO 


gusE 


jlOaj 


UELJ 


TOTHAO 


lOlQU 


QENV 


OTANA 


;o.o 


1 . 1 7 iE 02 


0.0 


-l.iflOE OA 


2.3A6E 05 


3.A5 IE 


OA 


9.26AE 06 


1.338E 05 


i.69WL 06 


i ODUNP 


kCjRK IN 


QABS 


ORE J 


HRS 


OOH 




OAUX 


TUTOL 




0.0 


3.399F C6 


•6.76<*E 06 


1.021E 07 


3.9O0E 02 


2.1336 


06 


1.032E 06 


-I.A20E 07 






Uni I 1 


WAS CALLED 


297B TIMES 
















2 




2970 — - 
















3 




9010 
















A 




6A0 7 
















5 




2978 
















A 




2976 
















7 




7A5 
















8 




32 














*" ♦ 


ENO OF RUN 


m 

















TiMf - 21:32250.62 
ELAPSE _T^IME_-_ lj.66 $EC. 
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Thes i s 

078877 

c.l 



Ounbar 



187979 



Analysis of a res- 
'dential heating 
system utilizing a 
solar assisted water- 
to-air heat pump. 
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Thes i s 
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Analysis of a res- 
•dential heating 
system utilizing a 
so^ar assisted water- 
o*r heat pump. 
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Analysis of a residential heating system 
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